Nonlinear spin-polarized transport through a ferromagnetic domain wall 
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A domain wall separating two oppositely magnetized regions in a ferromagnetic semiconductor 
exhibits, under appropriate conditions, strongly non linear I-V characteristics similar to those of a 
p-n diode. We study these characteristics as functions of wall width and temperature. As the width 
increases or the temperature decreases, direct tunneling between the majority spin bands reduces 
the effectiveness of the diode. This has important implications for the zero-field quenched resistance 
of magnetic semiconductors and for the design of a recently proposed spin transistor. 

PACS numbers: 



It has recently been reported that some doped semi- 
conductors, such as Gai^Mn^As (jj and Tii-zOo^O^ 
undergo ferromagnetic transitions at temperatures 
as high as 110 K and 300 K respectively, while others (n- 
doped Zni-^Mn^Se ^]) are almost completely spin po- 
larized by the application of a relatively modest magnetic 
field. These findings have raised hopes for the realization 
of semiconductor-based magnetoelectronic devices Q . 

In a ferromagnetic semiconductor, the up- and down- 
spin components of just one carrier type are quite analo- 
gous to majority and minority carriers in ordinary doped 
semiconductors. Accordingly, a domain wall separating 
two ferromagnetic regions with opposite magnetizations 
is the analogue of a p-n junction, while two consecutive 
domain walls correspond to a p-n-p transistor. In a re- 
cent paper || we have exploited this analogy to show that 
nonlinear amplification of a spin-polarized charge current 
is indeed possible in the "p-n-p" configuration, and can be 
controlled by a magnetic field or a voltage applied to the 
"base" region between the two domain walls. However, 
the analysis of Ref. || was based on the assumption that 
the probability of a carrier flipping its spin while crossing 
the domain wall is negligible. This corresponds to assum- 
ing the resistivity of the domain wall is large compared 
to that of the bulk material. 

The resistance of a domain wall between ferromagnetic 
materials has been examined several times from differ- 
ent perspectives since the pioneering work of Cabrera 
and Falicov Q. These authors found that the resistance 
was very small, and later calculations |7[ ||] have sup- 
ported that result for metallic magnets. A far differ- 
ent regime is possible, however, when the spin polariza- 
tion is or approaches 100%. For example, experimen- 
tal and theoretical results || indicate that domain walls 
in Lao.7Cao.3Mn03 may dominate the resistance in thin 
films. Magnetic semiconductor systems, due to their very 
small bandwidths, are also likely to be 100% spin polar- 
ized, and thus their domain walls should be highly resis- 
tive in the absence of spin-flip transport processes across 
them. 

A key question that has not been addressed so far is 



how the nonlinear current-voltage (I-V) characteristics of 
the domain wall are affected by spin-flip processes as the 
width of the domain wall increases. Note that the width 
of a domain wall can now be directly measured |ic|] and, 
in principle, geometrically controlled Our analytical 
theory of transport across the domain wall should there- 
fore be useful in designing devices with optimal values 
of the controllable parameters. Certainly such a theory 
would be crucial to understanding the zero-field quenched 
resistance and the low-field magnetoresistance of mag- 
netic semiconductors as well as to the realization of the 
"unipolar spin transistor" proposed in 

Here we present a quantitative study of the nonlin- 
ear I-V characteristics of a magnetic domain wall. The 
main issue is the competition between minority spin in- 
jection, which is responsible for the nonlinear spin-diode 
behavior, and majority spin transmission, which tends 
to suppress it. We shall show that the latter dominates 
when cither the temperature is low, or the domain wall is 
thick. Assuming that the motion of carriers through the 
domain wall is ballistic, we derive analytic expressions 
for the charge and spin currents as functions of applied 
voltage, width of the domain wall, and temperature. We 
further identify a new transport regime for intermediate 
wall thicknesses, in which carriers are ballistically trans- 
ported across the domain wall (characterized by nonlin- 
ear charge currents), but most spin polarization is lost. 

Our model is schematically depicted in Fig. 1(a). The 
two ferromagnetic regions FT and F2 are connected by 
a domain wall region of width d, —d/2 < x < d/2. The 
exchange field B(x) has the form 



B(x) = B [cos9(x)x + sm9(x)y], 



(1) 



where x, y are unit vectors in the direction of x and y, 
and the angle 9(x) varies linearly from 6 = tt/2 in F2 to 
6 = -n/2 in FT ||. 

We assume that d, while possibly large in comparison 
to a typical carrier wavelength, is smaller than the mean 
free path and the spin diffusion length L Sl which is in 
turn smaller than the geometric size of the system. A 
charge current J q is injected from the left: our objective 
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FIG. 1: (a) Schematic representation of a domain wall, (b) 
Qualitative behavior of the quasichemical potentials and the 
electrostatic potential (solid line). Note that the nonequilib- 
rium voltage drop occurs within the interfacial region, while 
the nonequilibrium populations extend up to a distance of or- 
der L s from it. (c) Reflection and transmission processes for 
an electron incident on the domain wall. 

is to calculate the voltage V that develops across the 
domain wall and the spin current J s due to the flow. 

Let /i> and /i< be the quasi-chemical potentials, which 
control the nonequilibrium densities of majority and mi- 
nority spin carriers respectively Far from the wall 
we have fi > = //< and the carrier densities have the equi- 



librium values n> and with riy » n 
variations from equilibrium A?i > ( < ) = «>(<) 
related to the difference of the quasichemical potentials 
A/i = /!< — fiy near the domain wall. Since, by charge 
neutrality, An < ~ — An> we see that the relative change 
in the minority spin density is always much larger than 
the corresponding relative change in the majority spin 
density. This implies that //> is essentially pinned to 
its bulk value, while /z< varies significantly in a region 
of length ~ L s on either side of the domain wall. We 
can therefore set /x> ~ throughout Fl and ~ eV 
throughout F2, where V is the electrostatic potential of 
Fl relative to F2 (see Fig. 1(b)) and the carriers are 
assumed to be electrons. The density variations are 
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where fcs is the Boltzmann constant and T is the tem- 
perature. 

The charge currents for majority and minority spin 
orientations must satisfy the condition J> + J< = J q 
where the total charge current J q is independent of po- 
sition. In addition, the minority carrier current J< is 
almost entirely a diffusion current, and is given by the 
classical relation J<{x) — eDdn < (x)/dx, where D is 
the diffusion constant. Because the spin density re- 
laxes to equilibrium exponentially on the scale of L s ( 



i.e., An < (x) — An < (±d/2)e~\ XTd / 2 ^ Ls where the lower 
sign holds in Fl and the upper sign in F2), the mi- 
nority carrier current at x — ±d/2 can be written as 
J<(±d/2) = TeDAn < (±d/2)/L s , or, with the help of 
Eq. (|), 



J<(±d/2) = T 



eDn 



(0) 



e An(±d/2)/k B T _ 1 



(3) 



It will be argued below that for nondegenerate carriers 
the quasi-chemical potential of minority spin electrons on 
each side of the domain wall adjusts to the quasichemical 
potential of majority spin electrons on the opposite side, 
so that j u < (-d/2) ~ eV, fi<(d/2) ~ (see Fig. 1(b)), 
and 



Afi(±d/2) = TeV. 



(4) 



Under the same assumption of nondegeneracy, it will 
also be shown that the matching condition for the spin 
current J s (x) = J-\{x) — J[{x) is 



J s (-d/2) _ t-+t+e 



-eV/k B T 



J s (d/2) t+ + t_e- cV / k BT 



(5) 



where t± — i n f ± i s f , and i s f and t n f are population- 
averaged transmission coefficients, with and without spin 
flip (see Fig. 1(c)), which will be defined more precisely 
below. Thus, the spin current is conserved across a sharp 
domain wall (t + = t_), but reverses its sign across a 
smooth one (i+ = — i—). 

Combining Eqs. (&p), and using current conservation 
we arrive at our main results. First 
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where Jq = 2eDn</L s . For t 



this reduces to 



the equation Q derived in (0), while for i n j = we 
get V = as expected for a ballistic conductor. In the 
linear regime eV/kgT << 1 this formula leads to the 
well-known interfacial resistance of Fert and Valet [ p5[ . 
Second, in the immediate vicinity of the domain wall the 
spin current is given by 



^ = 2sinh 2 (J*-) 
Jo \2k B Tj 



liktanhpM 



(7) 



where the upper sign holds in F2 and the lower sign in 
Fl. We see that spin-flip processes cause the appear- 
ance of an odd-in- voltage component of the spin-current, 
whereas, for t s f = 0, the spin-current is an even func- 
tion of F §. Shown in Fig. 2 is (a) the spin current 
in Fl, (b) the charge current, and (c) the ratio of the 
two. The curves correspond to several different values 
of t n f /t s f . The trends for the spin and charge current 
described above are evident in Fig. 2; specifically the 
charge current is always odd in V whereas the spin cur- 
rent is even in the absence of spin-flip. When spin-flip 
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FIG. 2: (a) Spin current in Fl, (b) charge current, (c) ratio 
of spin current to charge current vs. voltage for t n j /t s f = 
10, 2, 1, 0.5. 

dominates the spin current becomes odd as well. The 
spin current in F2 is related to that in Fl according to 
the following relation: J S (F2;V) = -J S (F1;-V). As 
tnf/isf becomes smaller, the "leakage current" between 
the two majority bands becomes significant, and the odd 
in V term in the spin current begins to dominate. Over 
the entire range shown of i n f /isf the relationship be- 
tween J q and V is highly nonlinear indicating ballistic 
transport. Thus ballistic transport itself is not a sufficient 
condition for maintaining spin polarization in transport 
across a domain wall. 

Assuming ballistic transport in the wall region, we cal- 
culate the transmission/reflection coefficients from the 
exact numerical solution of the Schrodinger equation 
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where A = gfiBB is the exchange spin-splitting. The 
technique of solution is the same as used in Ref. 0. 
Sample results are shown in Fig. 3(a)-(c) for three 
different values of the dimensionless parameter £ = 
Hn/2dy2mA = 10, 1, and 0.1, corresponding to sharp, 
intermediate, and smooth domain walls respectively, 
smooth domain walls respectively. Recent experiments 
p0| suggest the width of domain walls in artificial nanos- 
tructures can be as small as 1 nm, giving £ ~ 1 for an 
effective mass m equal to the electron mass and a spin 
splitting A = lOOmeV. Domain walls thinner than 20 nm 
have already been inferred in thin GaMnAs layers Eg- 

Fig. 3(d)-(f) shows the behavior of the key ratio 
tnf/isf as a function of temperature and thickness. As 
expected t n f vanishes at low temperature, because, in 
this limit, there are no incident states above the exchange 
barrier to provide minority spin-injection. The spin diode 
is a thermally-activated device (as a p-n diode is), thus 
higher temperature is favorable to its performance. Fig. 
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FIG. 3: (a)-(c) Energy dependence of transmission coefficients 
for £ = 10, 1, and 0.1 respectively . (d)-(e) Ratio of the pop- 
ulation averaged non-spin-flip to spin-flip transmission coeffi- 
cients {tnf/isf) vs. temperature for £ = 10 and 1 respectively, 
(f) W/V versus £ for k B T/A = 0.25 and 0.5. 

3(d,e) supports this view by showing that minority spin 
injection only dominates above a certain temperature 
(depending on domain wall thickness). However the con- 
dition ksT < A must be respected if the system is to 
be nearly 100% spin-polarized. The conclusion is that 
there is a range T m i n < T < T max in which unipolar spin 
diodes and transistors are expected to be operational. 

We now come to the justification of the matching con- 
dition (|^) and the calculation of the quasi-chemical po- 
tential offset. We begin with the former. In the spirit of 
the Landauer-Buttiker formalism we treat the ferromag- 
netic regions Fl and F2 as two reservoirs of spin polar- 
ized electrons at chemical potentials /ii = and [i 2 = eV 
which inject up- and down-spin electrons, respectively, 
in the domain wall region. The small density of minority 
spin carriers is neglected in the following argument. The 
components of the current due to electrons with energies 
in the range (E, E + dE) on the two sides of the domain 
wall are given (in units of e/h) by 

n> (E) = -(l~r nf (E))f 1> (E) + t sf (E)f 2> (E) 

31< {E) = r sf (E)f 1> (E) + t nf (E)f 2> (E) 

J2> (E) = (l-r nf (E))f 2> (E)-t. / (E)fi > (Ei) 

J2< (E) = -r sf (E)f 2> (E)-t nf (E)f 1> (E) 1 (9) 

where r ra / and r s f are the non spin-flip and spin-flip re- 
flection probabilities, related to t n f and t s f by the uni- 
tarity condition r n f + r sf + t nf + t s f = 1, and /i> ,/ 2 > 
are shorthands for the equilibrium distributions of ma- 
jority spin carriers in Fl and F2 respectively. Note 
that, for nondegenerate carriers /i> = f 2 >e~ eV / kBT . 
We find that the spin-flip reflection coefficient r s f is ex- 
tremely small at all energies and thicknesses, and can 



4 



therefore be safely neglected. With this approxima- 
tion, combined with the unitarity condition, it is easy 
to show that the energy-resolved currents are given by 
j sm (E) = (*_(+)(£) +t +{ _ ) (E)e~ eV / k * T )h > (E). Not- 
ing that f2>(E) oc e - E / k sT anc j integrating over energy 
we see that the total current J s i = / °° j s i(E)e 



-e/k B T 



IS 



equal to A(t_ + t + e eV / k sT^ wnere ^ ne average trans- 
mission coefficients are defined as 



t 
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J q °° e -E/k B T dE 



(10) 



and A is a constant. Similarly J s2 — A(t + +t-e~ eV / kBT ). 
The ratio J s \/ J S 2 is thus given by Eq. (g). 

To justify the quasi-chemical potential offset condi- 
tion, Eq. (^) we notice that the quasi-chemical potential 
/i<,i of minority spin electrons near the left hand side 
of the domain wall is an average of the quasi-chemical 
potentials of right (+) and left (-) moving electrons : 
e M<,i/fc B T = ^t^/ksT + e f*<,i/ fa}T ]/2. (A similar rela- 
tion holds for the quasi-chemical potential ^< ; 2 of minor- 
ity spin electrons near the right hand side of the domain 
wall). The quasi-chemical potentials for right and left 
movers on either side are determined by the conditions 
of continuity 



e -(B-/i+ 2 )/fc B T = 

e -(E-/i+ 2 )/fcsT = 

e -(E-t,- A )/k B T = 
e -(E-^- A )/k B T 



qe -(B-n+ A )/k B T + -(E-n% A )/k B T 
qe -{E-n+ A )/k B T + pe -{E-nt A )/k B T 



(11) 



where q = t sf / (t sf +t nf ) and p = t nf / (t sf +t nf ) are the 
relative probabilities of transmission with and without 
spin flip respectively. The first of these equations, for 
example, says that the density of right-moving up-spin 
electrons of energy E on the right hand side of the domain 
wall is equal to the density of right-moving down-spin 
electrons of the same energy which enter from the left 
and flip their spin, plus the density of right-moving up- 
spin electrons which enter from the left and do not flip 
their spin. Because the quasi-chemical potentials of the 
majority spin carriers are essentially pinned to their bulk 
values, we can set fit. 2 = fiZ 2 = H>,2 = eV and fit. ± = 
M> i = /- t >,i = 0- Integrating Eqs. ( p"l] ) over energy, and 
making use of p + q = 1, we easily get fi <t \ — eV and 
/i<,2 = 0, as indicated in Fig. 1(b). 

In summary, we have shown that both the thickness 
and the temperature have a profound influence on the 
nonlinear transport properties of a ferromagnetic domain 
wall. We have derived analytical formulas, Eqs. (^]) and 
(0), for the charge and spin currents of this "magnetic 
junction" under physical assumptions similar to the ones 
from which the Shockley equations of a classical p-n junc- 
tion are derived. These formulae indicate a new transport 
regime, where charge transport is ballistic, but spin po- 
larization is lost. Equations (||) and (0), together with 
microscopic calculation of the population-averaged trans- 
mission coefficients, can be used to assess the effective- 
ness of unipolar spin-diode devices in realistic circum- 
stances. 
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